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In a previous publication [1], an iterative closed-loop carrier synchronization 
scheme for binary phase-shift keyed (BPSK) modulation was proposed that was 
based on feeding back data decisions to the input of the loop, the purpose being 
to remove the modulation prior to carrier synchronization as opposed to the more 
conventional decision-feedback schemes that incorporate such feedback inside the 
loop. The idea there was that, with sufficient independence between the received 
data and the decisions on it that are fed back (as would occur in an error-correction 
coding environment with sufficient decoding delay), a “pure” tone in the presence 
of noise would ultimately be produced (after sufficient iteration and low enough 
error probability) and thus could be tracked without any squaring loss. This article 
demonstrates that, with some modification, the same idea of iterative information 
reduction through decision feedback can be applied to quadrature phase-shift keyed 
(QPSK) modulation, something that was mentioned in the previous publication but 
never pursued. 


I. Introduction 

In recent years there has been an ever-increasing interest in and application of highly power efficient 
error-correction codes such as turbo codes and low-density parity check (LDPC) codes to NASA-sponsored 
programs. These codes, which approach the Shannon channel capacity of the system, operate at very low 
symbol signal-to-noise ratios (SNRs), thus necessitating the need for synchronization (sync) (e.g., carrier, 
symbol, etc.) schemes that likewise operate efficiently at these SNRs. To begin to address this issue, 
several years ago the notion of iterative information-reduced carrier synchronization for coded binary 
phase-shift-keyed (BPSK) modulation was introduced [1], The idea behind this notion was to overcome 
the penalty in noisy reference loss attributed to the large squaring loss at low SNRs that is characteristic 
of the traditional types of BPSK carrier synchronization loops, e.g., the Costas or in-phase-quadrature 
phase (I-Q) loop, by removing the modulation from the received signal prior to demodulation, thus 
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allowing instead the use of a phase-locked loop (PLL) . The process by which the modulation was removed 
(information reduced) was envisioned as an iterative one wherein hard decisions on the data produced 
by the data detector/decoder would be fed back and applied to (multiplied by) the received signal, 
thus improving the carrier synchronization, which in turn would improve the data detector performance, 
hence the data decisions fed back on the next iteration, and so on. At the same time, the loop would be 
adapting its parameters (i.e., the loop nonlinearity in the in-phase arm) to match the current fidelity of 
the data decision process, all the while shifting its configuration from a Costas-type loop toward a PLL 
in the limit of perfect detection. Several publications based on this notion appeared in the literature 
and included everything from the basic idea and accompanying analysis/performance evaluations [1] to 
successful application and implementation for specific error-correcting codes [2,3]. 


With the continued shift toward increased bandwidth efficiency using modulations such as quadriphase- 
shift keying (QPSK), it became apparent, and was alluded to in the conclusions of [1], that the same 
ideas could be applied in this case, where the problem of nonlinear (e.g., quadrupling) loss (analogous 
to squaring loss for the BPSK case) at low SNRs becomes exacerbated. However, the manner in which 
the fed-back decisions are used for information removal is somewhat different since in addition to pro- 
ducing the desired (information-reduced) signal on, say, the I (or Q) channel, one now also generates an 
undesired (cross-talk) signal on the Q (or I) channel. A positive solution to this problem is proposed in 
this article and will be referred to as iterative QPSK information-reduced, carrier sync. As in the BPSK 
case, the scheme is motivated by maximum a posteriori (MAP) estimation of carrier phase considera- 
tions. Section II presents the derivation of the closed-loop architecture. Section III derives the tracking 
performance of the loop in terms of its mean-square phase error when operating in the linear (high loop 
SNR) region, as is typical. Section IV presents the asymptotic behavior of this performance in the limit 
of low SNR and compares it with that of the conventional QPSK tracking loop that does not make use 
of decision feedback for information reduction. Section V presents the numerical illustration of the loop 
tracking performance and compares it with that of the conventional loop mentioned above. Section VI 
explores the sensitivity of the scheme to mismatch between the actual value of symbol-error probability 
of the I and Q data decoders and the estimate of it used in the loop implementation. Finally, Section VII 
documents our conclusions. 


II. MAP Estimation of Carrier Phase for QPSK with Unbalanced Data 

Consider the transmission of a QPSK signal over an additive white Gaussian noise (AWGN) channel 
whereupon the received signal has the form x(t\ 9) = s(t; 9) + n(t), where 


s(t\ 9) = VS mi (t) sin(w c f + 9) + \/~Sm 2 (t) cos (uj c t + 9) (1) 

and the additive noise has the narrowband expansion 

n(t) = \/2[n c (f) cos ( u> c t + 9) — n s ( t ) sin ( uj c t + 0)] (2) 


In Eq. (1), S denotes the received signal power, 9 is random (assumed to be uniformly distributed and 
time independent) carrier phase, and mi (<) and m 2 (f) are the I and Q modulations, which are described 
by binary pulse trains as 
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oo 


m l(t) = ^2 CLkp(t — kT) 

k=— oo 

(3) 

OO 

m 2 (t) = bkp(t~kT) 

k=— oo 

where p{t) is a unit power rectangular pulse of duration T and {a*, },{&/,-} are each equiprobable inde- 
pendent, identically distributed (i.i.d.) binary (±1) sequences that are also independent of each other. 
The noise process in Eq. (2) is modeled in terms of a pair of independent, low-pass Gaussian processes, 
n c (t),n 8 (t), each with single-sided power spectral density (PSD) N 0 W/Hz and bandwidth Bh < w c /27t. 

With reference to the ad hoc structure illustrated in Fig. 1, the received signal-plus-noise x{t\ 9) is 
delayed by an amount A (to match the decoding delay of the I-channel data) and multiplied by the data 
estimator waveform 


OO 

rhi (t — A) = ^2 bkP ( t — kT — A) (4) 

k=— oo 

to form the signal yi(t) = x (t — A; 9) rhi (t — A). Similarly, the received signal plus noise x (t; 9) is shifted 
by — 7t/2 rad, delayed by an amount A (to match the decoding delay of the Q-channel data assumed to 
be the same as that of the I-channel decoding delay) and multiplied by the data estimator waveform 


OO 

rh 2 (t — A) = ^2 bkP ( t — kT — A) 

k=— oo 


( 5 ) 


to form the signal y 2 (t) = x (t — A; 9 — 7t/2) m 2 (t — A). Assuming, as in the previous publications on 
the subject, that the sequences of estimates {dk} , j^fc} are eac ^ i i d. sequences and have the probability 
statistics (based on error probability associated with the particular coding/decoding model employed) 


Pr {dk a-k} = Pi, Pr{a k =a k } = l~pi 


Pr jSfc yf 6 fc | = p 2 , Pr {& fc = & fc j = 1 - p 2 
then the results of the above-mentioned products become 

2 /i (t) = VSm i (■ t — A) mi (t — A) sin (u) c t + 9) + VSm 2 (t — A) rhi (t — A) cos ( uj c t + 9) + n\ (t) 

= VSe ii (■ t ) sin (u c t + 9) + VSe 2 1 (t) cos (tu c t + 9) + n\ ( t ) 

2/2 {t) = — VSmi (■ t — A) rh 2 (t — A) cos {w c t + 9) + VSm 2 (t — A) m 2 (t — A) sin ( oj c t + 9) + n 2 (t) 

= —VSe i 2 (t) cos (u c t + 9) + VSe 22 (t) sin (u c t + 9) + n 2 ( t ) 


( 6 ) 
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where we have arbitrarily absorbed the additional carrier phase, w c A, due to the decoding delay in the 
unknown carrier phase 0. The baseband modulation waveforms en (t) , e 22 (t ) , ei 2 (t ) , e 2 i (t) that reflect 
the errors in the detection of the I and Q input data sequences are modeled as 2 


where 


OO 

en (t) = ^2 e n kP(t-kT), e 2 2 (t) 

k =— oo 


OO 

ei 2 (t)= ^2 e i 2 kP(t-kT), e 2i (t) 

k =— oo 


oo 

e 2 2kP ( t - kT) 

k =— oo 


oo 

e 2 i fc p (i - fcT) 

k =— oo 


eiife = afchfe, e 22 fe = b^bk , ei 2 fc = afc&fc, e 2 ifc = cikbk 


(8) 


(9) 


are binary (±1) sequences with probability statistics 


Pr{e,:ife = -1} = K, Pr {e iifc = 1} = 1 - ft, * = 1,2 

Pr{ei 2 fc = -1} = Pr{ei 2fe = 1} = Pr{e 2U . = -1} = Pr{e 2 i*, = 1} = ^ 


(10) 


The additive noises in Eq. (7), namely, ni(t) = rh\(t — A)\/2 \n c (t) cos (ui c t + 0) — n a (t) sin(u; c f + 0)] and 
n 2 (t) = m 2 (t — A)v / 2[n c (t) sin (w c t + 0) + n s (t) cos {u) c t + 0)] are AWGN in each T-second interval with 
statistics identical to n(t) and are independent of the data. Note that to the extent that rh\ (t — A) and 
m 2 (t — A) are independent, n\ (t) and ?r 2 (t) are independent. 

In an effort to eliminate the quadrature-modulated carrier (cross-talk) term, we next form the sum of 
2 /i (t) and y 2 (t), resulting in 


y(t) = yi{t) + y 2 (t) 


= VS[ei(t) + e 2 (i)] sin(w c f + 0) + \Zs[e 2 i(t) - e 12 (t)] cos (w c t + 0) + ni{t) + n 2 (f) (11) 

= \fSEd{t) sin(w c f + 0) + VSE u (t) cos (w c t + 0) + ni(t) + n 2 (f) 

where the “d” subscript denotes the desired signal to be tracked and the “w” subscript denotes the residual 
undesired cross-talk signal; Edit) and E u (t) are ternary waveforms that can be expressed as 

OO 

Edit) = ^2 E dk.p{t - kT) 

k— — oo 

(12) 

OO 

E u {t) = ^2 E ukP{t ~ kT) 


2 Without loss in generality, for simplicity of notation we herein ignore the decoding delay. 
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where {E dk = o-kdk + bkbk} and {E uk = a k b k — a k b k } are ternary (0, ±2) i.i.d. sequences with probability 
distributions as follows: 


P{E dk } 


(1 -pi) (1 -p 2 ) , E dk = 2 


P 1 P 2 , = ^2 

Pi (1 -P 2 ) + P 2 (1 ~Pi) , E dk = 0 


(13) 


and 


( \vi (1 - P 2 ) + 7 .P 2 (1 - Pi) , E uk = ±2 
P{E uk }=l 2 2 (14) 

{ 1 -Pi (1 -P 2 ) -P 2 (1 -Pi) i E uk = 0 
For equal decoding probabilities, i.e., pi = P 2 = p, Eqs. (13) and (14) simplify to 


(1 -p) 2 , ^dfe = 2 

P{E dk }={p 2 , E dk = - 2 

2p(l-p), E dk = 0 


(15) 


and 


-P{E„fc} = 


’p(l -P) , 


E uk = ±2 


1 - 2p (1 - p) , E uk = 0 


(16) 


For p— 1/2, they have identical first-order statistics given by 


P {E dk } = P {E uk } = P {E k } = 


E k = ± 2 


-, E k — 0 


(17) 


Note that the sequences {E dk } and { E uk } are uncorrelated for all values of p but are not independent. 
Their joint statistics are given by 


P {E dkl E uk } 


' (1 - p) 2 , E dk = 2, E uk = 0 
P 2 , E dk = —2 ,E uk = 0 

p(l - p)i E dk = 0, E uk = ±2 
, 0, otherwise 


(18) 


Thus, when p = 1/2 (i.e., in the very early stages of iteration when the loop begins to operate without 
reliable data decisions), the signal component of the input to the loop as characterized in Eq. (11) appears 
as an I-Q modulation with identically distributed (but not independent) ternary modulations on each 
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arm. On the other hand, when p becomes small (the desirable case for reliable data detection), then in 
the limit as p — » 0, the signal component of the loop input becomes a pure (unmodulated) sinusoidal 
tone, which can be tracked by a PLL that exhibits the least degradation in performance. 

Returning now to the determination of the MAP estimate of phase, the combined signal plus noise in 
Eq. (11) is observed for K data intervals, i.e., over the interval 0 < t < KT. Based on this observation 
and knowledge of S,p(t) and uj c , the MAP estimate of phase is that value $map that maximizes the 
conditional probability density function p(9\y(t)) or, equivalently, the likelihood function p(y(t ) 1 0). As 
in previous analyses for equiprobable data statistics, the likelihood function (LF) conditioned on the 
data sequences {Edk} and {E uk } in the AT-symbol observation is given by (disregarding a proportionality 
constant) 3 


K- 1 

p(v (t) {Edk} , {E uk }) = U exp ( E dk y sk (9)) exp ( E uk y ck ( 9 )) exp 

k = 0 


E d ( 

y \^dk 


+ Eu k ) 


(19) 


where R d = ST /N 0 denotes the symbol signal-to-noise ratio (SNR) and 


a \fs r < ' k+1 '> T 

Vsk (0) = — / y ( t ) sin (u> c t + 9) dt 

1SI 0 JkT 


a %/S / > ( fc+1 ) T 

Vck (0) = - — y (t) cos (u c t + 9) dt 

iV 0 JkT 


(20) 


Jointly averaging Eq. (19) over the statistics of the data sequences { E dk } and {E uk } as given in Eq. (18) 
results in the unconditional LF 4 


k- 1 

A (6>) = p(y (t) \9 ) = ]^[ exp (-R d ) 

k—0 


(1 ~p) 2 exp (2 y sk (9)) + p 2 exp (-2 y sk (9)) 


+ 2p(l- p) cosh (2 y ck (6)) 


(21) 


Taking the natural logarithm of Eq. (21) gives the log-likelihood function (LLF) 


In A (6) 


k- 1 

- KR d + Y, ln 

k = 0 


(1 - p) 2 exp (2y ak ( 9 )) + p 2 exp (-2 y sk (9)) 


+ 2p(l- p) cosh (2 y ck ( 9 )) 


(22) 


The MAP phase estimate Omap is then the solution to 


3 Note that, because of the ternary nature of the effective I and Q data sequences {Edk} an d {E u k}i we temporarily retain 
the energy term in Eq. (19), exp [- (R d / 4) (Ej k + Ejfc)] > since, unlike the binary data case, it is possible that it is not 
independent of the data sequence. Shortly we shall show that indeed because of the dependence of Edk an d E u k on each 
other, the term E^ k + ^uk on iy t a k es on th e value 4, and thus this term may be dropped in the likelihood function. 

4 For simplicity, we pursue the equal decoding probability case. 
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(23) 


$map = max 1 In A(9) 

9 

Note that as the decoding error probability iteratively improves and eventually becomes vanishingly small, 
in the limit of p — > 0, Eq. (22) simplifies to (ignoring the constant additive term) 


A'-l 

In A (0) = 2 E Usk ( 0 ) 

k—0 


(24) 


which is the appropriate result for the LLF in the case of an unmodulated carrier. 

Since by definition $map is the value of 9 that maximizes the LLF, an equivalent statement is that 
(?map is the value of 9 at which the derivative of the LLF has zero value (and the second derivative is 
negative). For estimates 9 of 6 in the neighborhood of $map> the derivative of the LLF will be positive 
or negative in accordance with the sign of 9 — 0map; thus, this derivative can be used as an error signal 
in a closed-loop synchronizer to steer the loop in the direction of a locked condition corresponding to 
6 = 0map ■ Taking the derivative of Eq. (22) with respect to 9 gives 

— In A (<9) = 2 V 1 (1 - pf exp (2 y sk (6Q) - p 2 exp (-2 y sk (6Q) 

dd k=0 (1 - pf exp (2 y sk ( 9 )) + p 2 exp (-2 y sk ( 9 )) + 2p (1 - p) cosh (2 y ck ( 9 )) V ° k 


LTZ y sk (6 ) (25) 

+ 2p(l- p ) cosh (2 y ck (9)) 

which suggests a cross-over type closed loop with cross-coupled nonlinearities in its I and Q arms, where 
by “cross-coupled” is meant the fact that each one acts on both the I and Q amplified matched filter 
outputs. In particular, in accordance with Eq. (25), defining the two-dimensional nonlinearities 


K-l 


-*E 


2 p (1 - p) sinh (2 y ck 


k ~ 0 (! - p) exp (2 y sk ( 9 )) + p 2 exp (-2 y sk ( 9 )) 


and 


h (u, v) = 


(1 — p) 2 exp (2 u) — p 2 exp (—2m) 


(1 — p) exp (2m) + p 2 exp (—2 m) + 2p (1 — p) cosh (2m) 


h (w, v) = 


4p (1 — p) sinh (2v) 


(1 — p) exp (2u) + p 2 exp (—2 u) + 2p (1 — p) cosh (2v) 


(26) 


(27) 


then the closed-loop QPSK carrier synchronizer motivated by the MAP estimation procedure is as il- 
lustrated in Fig. 2, where the value of p to be used in forming the nonlinearities /i(m,m) and /2(m, v) is 
derived from an estimate p provided by the data detectors and is iteratively updated. We herein refer to 
this loop as the MAP estimation loop for QPSK with unbalanced data. Note that, for p — > 0, /i(m, v) — > 2 
and / 2 (m, v) = 0, in which case Eq. (25) simplifies to 


d9 


K - 1 

In A (0) = 2 E Uck (d) 

k—0 


(28) 


8 




9 


Fig. 2. Information-reduced carrier synchronization loop. 


which also can be immediately determined from Eq. (24). In principle, this is the same result that 
would occur for BPSK and suggests that, in the limit of vanishingly small decoding error probability, 
the closed loop behaves similarly to a phase-locked loop (PLL) since no quadrature arm is required in its 
implementation. 


Finally, for p = 1/2 where, as previously noted, the input to the loop appears as a modulation with 
uncorrelated but dependent I and Q ternary data streams with identical first-order statistics [see Eq. (17)], 
the nonlinearities of Eqs. (26) and (27) become 


fi (u, v ) 


2 sinh (2 u) 

cosh (2 u) + cosh (2v) 


h («, v) 


2 sinh (2v) — ft \ 

cosh (2 u) + cosh (2v) 1 


(29) 


which suggests a balanced (symmetrical nonlinearities) cross-over type of loop. 


III. Tracking Performance of the MAP Estimation Loop for QPSK with 
Unbalanced Data 

Consider the information-reduced MAP estimation loop for phase synchronization of QPSK with 
dependent unbalanced ternary I and Q data illustrated in Fig. 2. The analysis of the tracking performance 
of the loop in Fig. 2 begins by paralleling the approach taken in [4] for the MAP-motivated loop; however, 
there are some significant differences that arise fairly early in the development and inhibit completion of 
the analysis. First, the I and Q nonlinearities in each arm of the loop in [4], namely, hyperbolic tangent 
functions, are not cross-coupled, i.e., each depends only on the amplified matched filter output of its 
respective arm. Second, several of the steps in the analysis approach carried out in [4] depend upon the 
I and Q nonlinearities being an odd function of their argument, which is clearly not the situation here. 
Finally, the assumption of independent I and Q data sequences, as was the case in [4] where the input to 
the loop is truly the received QPSK signal, is another important difference between the analyses when 
carrying the statistical averages over these sequences needed to arrive at the parameters that characterize 
the loop’s performance. Despite these differences, we begin the analysis here in the same fashion as in [4] 
and proceed as far as one can go until the point is reached where certain simplifying assumptions must 
then be invoked. 

Assuming unit input I and Q phase detector (multiplier) gains and demodulation reference signals 
r c ( t ) = v^cos [ijj c t + 6) and r s (t) = \/2sin (u) c t + 0), then, after amplification by \JS/2/N 0 , the sample- 
and-hold outputs u k (</>) and v k {<f>) in the interval (k + 1) T <t< (k + 2)T are given by 


a \fsj2 f^ T ^Sj2 

uk{4>) = / y(t)r s (t)dt = 


N( 


0 JkT 


An 


T\j — ( E dk cos 4> - E uk sin 4>) - Ai sin <j> - N 2 cos <j> 


(30) 


Vk{4>) = 


\fS/2 


r(k+i)T 


A, 


0 JkT 


y{t)r c (t)dt = 


^Sj2 


An 


T\j — (E dk sin tj) + E uk cos <f>) - A 2 sin cj> + Ai cos 4> 
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where N x , N 2 are zero mean, independent Gaussian random variables 5 with variance a 2 Ni = cr^ = NqT = 
cr 2 and (j> = 6 — 9 is the loop phase error. Rewriting Eq. (30) in normalized form, we have 


u k {4>) 


Vk ( 0 ) 


1 

2 

1 

2 


(-Bdfc cos (/> - E uk sin <j>) - \j2R d X x sin (j> - sj2R d X 2 cos </> 
Rd ( E dk sin (/> + E uk cos </>) - \/2R d X 2 sin + ^2R d X 1 cos 


(31) 


where X., = Ni/<r,i = 1,2 are normalized (unit variance) i.i.d. Gaussian random variables (RVs). Mul- 
tiplying v k ((f>) by the nonlinear ly processed u k ((f>),v k ((f>) pair gives one component of the dynamic error 
signal, namely, 


z k i{4>) = v k {4>)fi{u k (c/>) ,v k (<j>)) (32) 

Similarly, multiplying u k (4>) by the nonlinearly processed u k (</)) , v k (</>) pair gives the other component 
of the dynamic error signal, namely, 

z k 2 {<t>) = u k ((/>) f 2 (u k (<f>) ,v k (<j>)) (33) 

Finally, the total dynamic error signal, e k (</>), in the fcth signal interval is the difference of z k i((t>) and 
z k 2 { 4 >), i-e., e k ((f>) = z k \{<t>) - ^ 2 (</>)• 

The tracking performance of a loop such as that in Fig. 2 can, in its linear region of operation (small 
phase error), be characterized by the variance of the phase error which can be evaluated from the slope 
of the equivalent S-curve at <j> = 0 and the power spectral density (PSD) of the equivalent additive 
noise. These parameters in turn can be determined by evaluating the statistics of the signal and noise 
components of z Q (t; (f>), making the usual assumption that the loop bandwidth is much less than the data 
bandwidth. With this in mind, the loop S-curve is defined as the expected value of the error signal and 
is given by 


V O) = e k (< !> ) 


■Edk,E uk ,X i,X2 


= v k (</>) fl (lift (</>) , V k (4>)) - u k (</>) f 2 (u k {(/)) , V k (0)) 


■Edk, Euk ->X 1,X2 


(34) 


where the overbar denotes statistical averaging over the probability distributions of the data and noise RVs 
indicated, keeping in mind that the average over the data RVs must be performed jointly in accordance 
with Eq. (18). While, in principle, the averaging required in Eq. (34) cannot be obtained in closed form, 
using the properties of the nonlinearities obtained from Eqs. (26) and (27), namely that 


fi {u, v ) = fi (■ u , - v ) 
f 2 (u,v) = -f 2 (u,-v) 


(35) 


it is nevertheless possible to show analytically that the S-curve is an odd function through the ori- 
gin, i.e., 77 ( 0 ) = — 77 (— 0), as is desirable of such a tracking characteristic. Figures 3(a) through 3(c) are 


5 Note that multiplication of a Gaussian process by a binary (d=l) waveform of rate 1/T does not change the Gaussian 
properties of the original process within any T-second time interval. Thus, for any input bit, the corresponding sample- 
and-hold outputs of the integrate- and-dump (I&D) filters are still Gaussian random variables. 
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illustrations of the transitional behavior of the S-curves for <j> > 0 (obtained by computer evaluation 
of the averages required in Eq. (34) using MATLAB™ software) for three different values of p with 
symbol SNR Rd as a parameter. The y-axis scale on the left applies to the curves corresponding to 
the lower values of SNR (the dashed curves), whereas the y-axis scale on the right applies to the curves 
corresponding to the higher values of SNR (the solid curves) . In addition, we remind the reader that in 
the limit of p — > 0 the loop becomes a PLL, and thus the S-curve has the form of sin </>. 

The slope K v of this S-curve at the origin is obtained from 

K v = U=° = u k (0)fi(u k (0),Vk(0)) +Vk(0)f[(uk(<j>),Vk(<j>)) 1 4,=° 


E dk iE u k ,Xi ,-X^2 

+ v k (0)f 2 (uk(0),Vk(0)) - Uk(0)f!i(uk(</>),Vk(</>)) U =0 (36) 


where from Eq. (31) we have made use of the relations 


du k (</>) 
d(p 


= - v k (</>) 


dvk (4>) 

d<j) 


= u k (</>) 


(37) 


Also, the derivatives of the nonlinearities as denoted by the primes in Eq. (36) are given by 


f[ (u k (<£) , v k {&)) = 


dfi (u k (</>) , v k ((/>)) du k (0) dfo (u k (</>) , v k ((/>)) dv k (<j)) 


du k ((f)) 


d(j> 


dv k (</>) 


d(j> 


= ~Vk{<t>) ,,, + U k (<j>) a! 7T\ . t = 1 ’ 2 


du k (</>) 


dv k (</)) 


(38) 


From Eqs. (26) and (27), the partial derivatives needed in Eq. (38) are given by 


a dfi (u, v) df 2 (u, v) 
ff(u,v)= - 


du 


dv 


8p(l — p) |2p(l — p) + cosh (2v) (1 — p) 2 exp (2u) + p 2 exp (— 2u) | 


(1 — p) exp (2 u) + p 2 exp (—2 u) + 2p (1 — p) cosh (2v) 


, , a dfi (u, v ) (9/ 2 (u, v) 

h{u, v) = — 


dv 


du 


(39) 


8 p (1 — p) sinh (2v) (1 — p) exp (2 u) — p 2 exp (—2 u) 


(1 — p) 2 exp (2m) + p 2 exp (—2m) + 2p (1 — p) cosh (2m) 
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which for p = 1/2 simplify to 


g ( u , u) = 


dfi {u, v ) 5/2 (m, v) 


dii 


dv 


4 [l + cosh (2 m) cosh (2 k) 
[cosh (2m) + cosh (2m) 2 


^ ^ a d/i (it, it) 3/2 (m, m) 


4sinh (2m) sinh (2m) 
[cosh (2m) + cosh (2m) 2 


(40) 


and also could be determined directly by differentiation of Eq. (29). Finally, substituting Eq. (38) 
combined with Eq. (39) into Eq. (36), the slope of the S-curve becomes 


- - Z~Edk i^uk )-X”l 1X2 

K v = u k (0) /1 (u fe (0) , v k (0)) + M fc (0) f 2 (u k (0) , v k (0)) 

~r o ~/ \ O / 7 ' ", 7 . \ Edk T^uk j-^"l ,x 2 

- [u k (0) + v 2 k (0)]g(u k (0) , v k (0)) 

Z~Edk ,E uk 1X2 

+ 2 Mfc (0) Mfc (0) h(uk ( 0 ) , Mfc (0)) 

where, from Eq. (31), 


u fc(0) = 2 ( ' R-dEdk - V2RdX 2 ) 
v k ( 0 ) = - {RdE u k + \/2RdXi ^ 


(41) 


(42) 


We note that, in an actual implementation, p would itself depend on the loop phase error (since the data 
decisions being fed back to the input of the loop depend on the phase of the in-phase carrier demodulation 
reference derived from the loop). However, since the slope of the S-curve is evaluated at <f> = 0, then the 
value of p needed in Eq. (41) is the value at (f> = 0, i.e. , the ideal performance of the data detector. Thus, 
the dependence of p on (j> is irrelevant to the evaluation of the S-curve slope. 

The noise component of e k {t\(j>) (evaluated at <j> = 0) is given by N e (t) = e k (t; 0) — 77 ( 0 ) = e k (t\0), 
which from Eqs. (32) and (33) becomes 


N e (t) = v k (0)/i (ufc (0) ,M fc (0)) - Mfc (0) f 2 (u k (0) ,M fc (0)) (43) 


with variance 


~Edk ,E u k ,Xi ,X2 

= [ v k (0) /l(Mfc (0) ,Mfc (0)) - Mfc (0) / 2 (w fc (0) ,Mfc (0))] (44) 

Here again the value of p needed in Eq. (44) is the value at 4> = 0, which is consistent with the previous 
usage of the same parameter for the S-curve slope in Eq. (41). 

Because of the sample-and-hold circuits shown in Fig. 2, the noise process of Eq. (43) is piecewise 
constant over intervals of T-second duration. Thus, as long as the loop bandwidth is much less than the 
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data bandwidth, this process can be approximated, as has been done in the past, by a delta-correlated 
process with correlation function given by 


Rn s (t) = N e (t)N e (t+r ) 




M < T 


o, \t\>T 


and equivalent single-sided noise spectral density (in the neighborhood of / = 0) 


(45) 


/ OO 

R Ne (r)dT = 2a 2 N T 

-OO 


(46) 


Finally, then, the linearized phase error variance is given by 


a 


2 

4 > 


KBl 

K 2 


A 


(pSl) 1 


(47) 


where p = S/NqBl is the linear loop (PLL) signal-to-noise ratio (Bl is the single-sided loop bandwidth) 
and, analogous to conventional quadriphase Costas loop terminology, Sl is the “quadrupling loss,” which 
reflects the penalty paid due to the signal and noise cross-products present in Solving for Sl 

from Eq. (47) gives 


S L 


K 2 


Rd(N>/T) 


(48) 


which unfortunately cannot be obtained in closed form for arbitrary values of p. In the limiting case of 
p = 0, however, K v = 2 Rd and N' 0 = 4 RdT, and thus Sl = 1, which corresponds to no quadrupling loss, 
as one would expect of a PLL. 


IV. Asymptotic Behavior at Low SNR 

It is of interest to examine the behavior of the information-reduced carrier synchronization scheme at 
low SNR (small symbol energy-to-noise ratio) since this is the region where the quadrupling loss associated 
with conventional carrier synchronization schemes is large and thus limits their performances. Since the 
outputs of the I&D filters in Fig. 2 are proportional to Rd, then for small values of this parameter, the 
I and Q nonlinearities can be approximated by the first few terms of their Taylor series expansion. In 
particular, for the p = 1/2 case, making the approximations sinli x = x + x 3 /6, cosh a; = 1 + x 2 /2 in 
Eq. (29), we obtain (keeping only terms up to order a; 3 ) 


fi{u,v) = 2 



f 2 (u,v) = 2 



which when substituted in Eq. (34) results in the S-curve 


(49) 
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A E dk >^uk >-^1 > ^2 A R'dk ’^uk >-^1 >-^2 

V (<t>) = 3 K- (</>) *4 (</>) - Wfe (</’) <4 (</>)] = 3 Wfc (</>) Vk { 4 >) [tifc (</>) - V 2 k (</>)] (50) 

It is interesting (although not surprising) that the functional form of the error signal in Eq. (50) is identical 
to that in a conventional QPSK Costas loop (derived from MAP phase estimate considerations at low 
SNR) [4]. Substituting Eq. (31) into Eq. (50) and carrying out the required statistical averages results in 

V (</>) = \Ri sin 4^ (51) 

with slope at the origin 

K v = (52) 

To evaluate the equivalent noise power spectral density in Eq. (46), we need to evaluate the variance 
in Eq. (44), which, using the approximate expressions for the nonlinearities in Eq. (49), becomes 



However, from the definitions of Uk (0) and Vk (0) in Eq. (31), these variables are noise-wise independent. 
Thus, the following first-order moments are useful in evaluating the joint moments required in Eq. (53): 


tfM X2 = 1 -[R 2 d E 2 dk + 2R d ] 
vf{0) Xl = l -[R 2 d E 2 uk +2R d ] 

<W X ' 2 = ^ [RdE A dk + 12 R 3 d E 2 dk + 12R 2 ] 

(54) 

4¥) Xl = ^ [RiKk + 12 B? d Elk + 12 R 2 d ] 

<(0) X2 = l~ A [R G A + 30 R 5 d Ej k + 180 RA + 120 <R 3 d \ 

^0) Al = 1 [RA + SOi^ + 180 R\E 2 uk + 120 R 3 d \ 

Using Eq. (54), the joint averages over the data sequences of the terms in Eq. (53) become 
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( 55 ) 


4 (oK(o) 


E dk’ E uk’ X 1> X 2 


= «2(oK(o) 


E dk’ E uk’ X l’ X 2 


= jg (4-Rd + 30i$ + 607?* + 15 R 4 ) 


v 4 k (0) u 4 (0)^— 2 = ^ (I 2 i?e + 36 R 5 d + 9R 4 ) 


which upon substitution in Eq. (53) gives 

= l (8 Rl + 36 R 6 d + 48 R 5 d + 12 R 4 ) (56) 

Finally, using Eq. (56) in Eq. (46) together with Eq. (52), the quadrupling loss as computed from Eq. (48) 
becomes 


S L = 


1 

' 9 6 3“ 

+ 2R d + R^ d + 2R^ d 


(57) 


which is identical to the result for the conventional QPSK Costas loop derived from MAP carrier phase 
estimation considerations at low SNR [4, Eq. (35)]. What is particularly interesting about this result is 
that, despite the fact that for p = 1/2 the effective I and Q sequences that are input to the iterative 
information-reduced loop are ternary in nature and not independent of each other, whereas in the con- 
ventional QPSK Costas loop they are i.i.cl. binary and independent, the quadrupling loss is nevertheless 
the same in the limit of low SNR operation of the former. 

Thus, we see that at the outset of the iteration process, the information-reduced QPSK loop performs 
the same as the conventional QPSK loop. However, as the loop continues to produce better and better 
carrier phase estimates, p continues to decrease, and thus the quadrupling loss performance of the loop 
varies between its value as given by Eq. (57) and its value corresponding to the error rate performance 
of the data detector at the particular value of R d under consideration. (In the limit of p — > 0, the 
quadrupling loss would disappear entirely.) By contrast, the quadrupling loss of the conventional QPSK 
loop remains fixed at its value as determined from Eq. (57). In the following section, we will demonstrate 
this more quantitatively with numerical results obtained by computer simulation. 


V. Evaluation of the Tracking Performance as a Function of Detection SNR 

As discussed in the previous section, the amount of quadrupling loss is commonly used as a measure of 
the tracking performance degradation of QPSK carrier synchronization loops. With this in mind, Fig. 4 
is a plot of the quadrupling loss Sl of the information-reduced QPSK loop versus detection (symbol) 
SNR R d with p as a parameter. These results were obtained by computer simulation (MATLAB™ 
numerical methods) of the statistical averages required to compute the S-curve slope K r) and equivalent 
noise PSD Nq needed in Eq. (48). While we recognize that, even for perfect carrier synchronization, 
any given form of coded modulation and associated data detector/decoder results in a value of p that 
is a function of R d , we have chosen to maintain p as a constant parameter in this figure. The reason 
for this is that the results presented there can then be independent of the specific form of data detector 
or, equivalently, the specific form of error correction coding employed, provided that the conditions for 
an i.i.d. and noise-dependent error sequence are still met. For any given coding application, only a sin- 
gle point on each curve of constant p would be applicable, namely, the one corresponding to the given error 
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Fig. 4. Quadrupling loss versus symbol SNR in decibels 
with p as a parameter. 


rate behavior of the data detector (decoder), and thus the performance for that particular application 
would be illustrated by a single curve passing through this series of points. 

We observe from the results in Fig. 4 the monotonic behavior of the quadrupling loss as a function 
of p. 6 Although somewhat cumbersome to evaluate analytically, the limiting behavior of this loss at low 
SNR is given by 


Jim S L = (1 - 2pf (58) 

rt-d — ►U 

which is the identical result obtained for the information-reduced BPSK loop in [1], Furthermore, we 
observe the immediate large improvement in quadrupling loss at low SNR as soon as some degree of 
fidelity of the data decisions begins to take place. More specifically, for improvement to take place in 
the carrier synchronizer, the data estimator is not required to operate with symbol-error probabilities as 
small as those typically needed for reliable communication. In fact, when used in the manner illustrated 
in Fig. 1, a value of p = 0.1, which in most applications would be considered unacceptable for data 
decisions, can result in dramatic improvement in carrier synchronization performance relative to that of 
the conventional QPSK Costas loop as represented by the curve labeled p = 0.5. It should be further 
noted that, although this curve was obtained by computer simulation, it agrees with the theoretical result 
in Eq. (57). 

Another interpretation of the numerical results can be obtained by recognizing that the loop SNR 
p = S/N 0 Bp can be expressed in terms of the symbol SNR Rd = ST /No by p = Rd/BpT , which when 
substituted in Eq. (47) gives an expression for the normalized phase error variance as 

2 

= ( R d S L r 1 (59) 


Figure 5 is a plot of this normalized variance, with Sl determined as in Fig. 4 versus Rd with p as 
a parameter. As one would anticipate, here the curves are both monotonic with p for fixed Rd and 


6 At the same time, we note that for fixed p we cannot guarantee a monotonic performance with SNR. 
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Fig. 5. Normalized phase error variance versus symbol 
SNR in decibels with p as a parameter. 

monotonic with Rd for fixed p. These curves represent the tracking performance of the loop when the 
loop bandwidth-symbol time product is held fixed and once again show the rapid dramatic improvement 
with data detection efficiency at low SNR for the information-reduced QPSK loop. Here again we see the 
dramatic improvement in performance at low SNR as the information reduction increases, i.e. , the data 
estimates improve. 

VI. Sensitivity to Mismatch 

Having improved the phase tracking by the use of data feedback, the new improved carrier phase 
estimate will help the data detector obtain better symbol estimates, which suggests an iterative procedure. 
That is, the phase tracking loop is initially designed assuming a value of p corresponding to no feedback. 
The initial data estimates obtained using this particular structure as a carrier synchronizer are fed back, 
and the carrier synchronizer structure then is modified based on the new estimated value of p. This 
updated carrier synchronizer provides improved phase tracking, which in turn provides improved data 
detection. These new data estimates, having a smaller value of p, then are fed back to the input of the 
carrier synchronizer and the iteration continues. Thus, in practice, at any given time instant there will be 
a mismatch between the true value of p associated with the input error sequence and its estimate, p , used 
to implement the nonlinearities of Eqs. (26) and (27) in the loop. As such, we wish to investigate the 
sensitivity of the loop performance as measured by the quadrupling loss to mismatch between p and p. 
That is, how good must our estimate of p be for the purpose of loop implementation and yet still achieve 
an improvement in performance relative to not implementing decision feedback at all? 

We begin this discussion by recognizing that when p and p are unequal, the quadrupling loss expres- 
sion of Eq. (48) is still appropriate provided that, in the expressions used to compute K v and Nq, p is 
substituted for p in the definitions of the nonlinearities in Eqs. (26) and (27) as well as their derivatives 
in Eq. (39). The first point to observe is that, for p = 0.5, f±(u,v) and f 2 (u,v) and their associated 
derivatives simplify to the forms in Eqs. (29) and (40), respectively, which in view of the symmetry prop- 
erties of these functions results in the averages over the data sequences E u p and Edk needed to evaluate 
K v and Nq becoming independent of p. This in turn results in the quadrupling loss also becoming in- 
dependent of p. That is, if we assume p = 0.5 for the loop design, then the performance is independent 
of the true data input probabilities. Next, suppose that we are optimistic in our estimate of p used for 
the design of the loop, i.e., p < p. Figures 6(a) and 6(b) are plots of the squaring loss as determined 
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Fig. 6. Quadrupling loss performance in the presence of 
mismatch: (a ) p = 0.0\,p<p and (b ) p = 0.1, p <p. 


from Eq. (48) versus Rd in decibels with p as a parameter and p = 0.01 and p = 0.1, respectively. The 
deviation between p and p is as much as a factor of 5. We see that for this case the performance of the 
loop is quite insensitive (a variation on the order of 0.1 dB or less) to the mismatch and yet we still 
obtain a significant improvement relative to the case of no feedback at all. Figures 7(a) and 7(b) show 
the analogous results for the case where we are pessimistic in our estimate of p used for the design of 
the loop, i.e. , p > p. Although we see a bit more sensitivity of the loop performance to the mismatch 
between p and p , in all cases performance is as good or better than the case of no feedback (information 
reduction) at all. It is also interesting to note in these figures that, except for the singular case of p = 0.5, 
the limiting quadrupling performance as Rd — > 0 behaves like (1 — 2 p) independently of the amount of 
mismatch and thus is consistent with Eq. (58) for the perfectly matched case. For p = 0.5, as previously 
mentioned, the performance is independent of the true value of p and becomes synonymous with that of 
the conventional Costas loop. 
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Fig. 7. Quadrupling loss performance in the presence of 
mismatch: (a ) p = 0.01, p^p and (b) p = 0A,p>p. 


VII. Conclusions 

Using the same iterative information-reduction principle previously introduced for carrier synchroniza- 
tion of coded BPSK, it is possible to design an analogous closed-loop configuration for coded QPSK that 
outperforms the conventional I-Q schemes. It has been shown that very accurate symbol decisions are 
not required in order to achieve significant improvements in low SNR tracking threshold (or, equivalently, 
operation with greatly reduced phase error) as measured by the reduction in the quadrupling loss of the 
loop. Stated another way, the new QPSK carrier synchronization loop can operate effectively at signal 
levels that might be too low for a conventional I-Q loop to even lock. Furthermore, once phase-locked, 
the loop continues to iteratively improve its performance by transferring its phase estimates to the data 
detector and receiving in return improved data decisions to modify its I and Q nonlinearities, all the while 
attempting to approach the tracking performance of a phase-locked loop operating on an unmodulated 
tone. 
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Finally, it should be noted that the same information-reduction principle can in principle be applied 
to offset QPSK (OQPSK) to generate a carrier tracking loop for such a modulation. However, in this 
case the effective ternary I and Q sequences generated at the input to the loop now occur at twice 
the data rate (because of the half-symbol offset between the transmitted I and Q data streams), and 
furthermore these effective data sequences are now themselves correlated in time in addition to being 
jointly dependent. Because of this, analysis of the behavior of such a loop is difficult if not impossible, 
and thus the performance would need to be obtained by computer simulation. Such a task is left as an 
exercise for the interested reader. 


References 

[1] M. K. Simon and V. A. Vilnrotter, “Iterative Information-Reduced Carrier Syn- 
chronization Using Decision Feedback for Low SNR Applications,” The Telecom- 
munications and Data Acquisition Progress Report 42-130, April-June 199 7, Jet 
Propulsion Laboratory, Pasadena, California, pp. 1-21, August 15, 1997. 

http: //ipnpr/progress_report / 42-130 /130A.pdf 

[2] V. Vilnrotter, A. Gray, and C. Lee, “Carrier Synchronization for Low Signal-to- 
Noise Ratio Binary Phase-Shift-Keyed Modulation Signals,” The Telecommuni- 
cations and Mission Operations Progress Report 42-139, July-September 1999, 
Jet Propulsion Laboratory, Pasadena, California, pp. 1-16, November 15, 1999. 
http://ipnpr/progress_report/42-139/139I.pdf 

[3] V. Vilnrotter, C. Lee, and N. Lay, “A Generalized Pre-Processor for Block and 
Convolutionally Coded Signals,” The Telecommunications and Mission Opera- 
tions Progress Report. 42-144 < October December 2000 , Jet Propulsion Labora- 
tory, Pasadena, California, pp. 1-19, February 15, 2001. 

http: //ipnpr/progress_report / 42-144/144G.pdf 

[4] M. K. Simon, “On the Optimality of the MAP Estimation Loop for Carrier Phase 
Tracking BPSK and QPSK Signals,” IEEE Transactions on Communications, 
vol. COM-27, no. 1,, pp. 158-165, January 1979. 


22 



